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(57) Abstract 

A direction for a desired radio signal 
arriving at a receiving station is determined 
from a plurality of estimated signal medium 
response coefficients for a ray of the desired 
radio signal. According to an aspect of the 
present invention, the estimated signal medium 
response coefficients for the ray are generated 
according to an interference rejection combin- 
ing (IRC) process. According to another as- 
pect, a direction of arrival for a ray of the de- 
sired signal is determined and used to compute 
a power spectral density for the ray for a plural- 
ity of directions. A direction of arrival for the 
ray is then determined by determining a direc- 
tion for which the power spectral density has a 
maximum value. In addition to desired signal 
direction, direction of arrival for an interfering 
signal can be determined using an impairment 
autocorrelation matrix that is also generated as 
part of the IRC process. 
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APPARATUS AND METHODS FOR DETERMINING SIGNAL DIRECTION 
FROM AN ESTIMATED SIGNAL MEDIUM RESPONSE FOR A RAY 
COMPONENT OF A RADIO SIGNAL 

Field of the Invention 

The present invention relates to communications apparatus and methods, more 
particularly, to apparatus and methods for determining a direction of a radio signal. 

5 Background of the Invention 

Wireless communications systems such as cellular radiotelephone systems typically 
are subject to performance-degrading environmental effects such as multi-path fading 
and interference. Fading effects include flat fading, arising from the interaction of a 
transmitted signal (the main ray) with reflected versions of the transmitted signal that 

1 0 arrive concurrently at a receiver. Time dispersion, another type of fading, arises from 
interaction of the main ray with time-delayed reflections of the main ray. Interference 
effects may be caused by interaction of non-orthogonal signals generated in the signal 
medium by sources other than the source of the desired transmitted signal. Well- 
known techniques for reducing the effects of flat fading include diversity combining 

1 5 of signals from spatially separated antennas. Equalization techniques such as 

maximum likelihood sequence estimation (MLSE) may be used to compensate for 
time dispersion. Interference may be reduced by using antenna beam steering to 
reduce reception of imdesired signals. 

United States Patent Application Serial No. 08/284,775, filed August 2, 1994, 

20 assigned to the assignee of the present invention and incorporated herein by reference, 
has proposed an interference rejection combining (IRC) technique which compensates 
for flat fading, time dispersion and interference. An impairment autocorrelation 
matrix is computed and used to compute branch metrics for a sequence estimation 
process which are corrected for the presence of interfering signals, i.e., impairments, 

25 using the autocorrelation properties of the interfering signals. According to United 
States Patent Application Serial No. 08/655,930, filed May 31, 1996, assigned to the 
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assignee of the present application and incorporated herein by reference in its entirety, 
IRC may be used to determine a covariance matrix of a desired signal, and downlink 
beamsteering techniques may be employed based on the determined covariance to 
improve communications between a mobile unit and a base station 
5 Unfortunately, however, estimates of direction of arrival of a desired signal 

produced using the approach described in the aforementioned United States Patent 
Application Serial No. 08/655,930 can be inaccurate due to noise introduced in 
determining the covariance matrix of the desired signal. The noise susceptibility may 
in turn require increased sampling to achieve a desired accuracy for the estimate of 
10 the direction of arrival. Accordingly, there is a need for improved techniques of 

direction estimation that are less susceptible to noise and can provide more accurate 
results without requiring undue sampling. 

Summary of the Invention 

15 In light of the foregoing, it is an object of the present invention to provide 

apparatus and methods for determining a direction for a desired signal that are 
resistant to fading and interference effects. 

It is another object of the present invention to provide apparatus and methods 
for determining a direction for a desired signal that can provide greater accuracy than 
20 conventional techniques. 

It is yet another object of the present invention to provide apparatus and 
methods for determining direction of an interfering signal. 

These and other objects, features and advantages are provided according to the 
present invention by apparatus and methods in which a direction for a desired radio 
25 signal arriving at a receiving station is determined from a plurality of estimated signal 
medium response coefficients for a ray of the desired radio signal. According to an 
aspect of the present invention, the estimated signal medium response coefficients for 
the ray are generated according to an interference rejection combining (IRC) process, 
an equalization process that compensates for interference. According to another 
30 aspect, a direction of arrival for a ray of the desired signal is determined and used to 
compute a power spectral density for the ray for a plurality of directions. A direction 
of arrival for the ray is then determined by determining a direction for which the 
power spectral density has a maximum value. By determining desired signal direction 
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Utilizing coefficients for a ray which are normally generated as part of an IRC 
process, a technique for determining signal direction is provided which is potentially 
more efficient than conventional techniques. In addition to desired signal direction, 
direction of arrival for an interfering signal can be efficiently determined using an 
5 impairment autocorrelation matrix that is also generated as part of the IRC process. 

In particular, according to the present invention, a direction for a 
desired radio signal generated in a communications medium is determined. A 
plurality of radio signals is received from the communications medium, e.g., at a like 
plurality of antennas. A plurality of estimated signal medium response coefficients for 
10 a ray of the desired radio signal are generated firom the plurality of received radio 
signals. A direction for the desired radio signal is determined fi-om the plurality of 
estimated signal mediimi response coefficients for the ray of the desired radio signal. 
A direction for the desired radio signal may be determined by determining a direction 
of arrival for the ray from the corresponding pluraUty of estimated signal medium 
1 5 response coefficients for the ray, and then determining a direction for the desired 
radio signal firom the determined direction of arrival for the ray. 

According to an aspect of the present invention, the plurality of estimated 
signal medium response coeflficients for the ray may be generated by generating a first 
plurality of estimated signal medium response coefficients. An impairment 
20 autocorrelation matrix is then generated from the plurality of received signals and the 
first plurality of estimated signal medium response coefficients. A branch metric is 
generated &"om the generated impairment autocorrelation matrix and the first plurality 
of estimated signal medium response coefficients, and a symbol estimate is generated 
firom the generated branch metric. A second plurality of estimated signal medium 
25 response coefficients may then be generated from the plurality of received signals and 
the generated symbol estimate. 

According to another aspect of the present invention, the plurality of estimated 
signal medium response coefficients may be generated by generating a series of 
pluralities of estimated signal medium response coefficients over a time interval, 
30 processing the series of pluralities of estimated signal medium response coefficients 
with a tracking filter to produce a plurality of predicted signal medium response 
coefficients, and determining a determining a direction for the desired radio signal 
fi"om the plurality of predicted signal medium response coefficients for the ray of the 
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desired radio signal. The desired radio signal may be formatted according to a slotted 
communications protocol defined according to a plurality of time slots, and the series 
of pluralities of estimated signal medium response coefficients may be generated over 
at least one of the plurality of time slots. 
5 According to another tracking aspect, a plurality of impairment autocorrelation 

matrices may be estimated over a time interval, e.g., over a time slot or a plurality of 
time slots in a slotted communications system. The generated plurality of impairment 
autocorrelation matrices may then be processed with a tracking filter to produce a 
predicted impairment autocorrelation matrix. A branch metric may then be generated 
10 from the predicted impairment autocorrelation matrix, and used to generate an 
improved direction estimate. 

According to yet another method aspect of the present invention, a direction of 
arrival for the ray is determined by determining an autocorrelation value for the ray of 
the desired radio signal. A plurality of power spectral density values may be then be 
1 5 determined for the ray at a like plurality of directions using the determined 

autocorrelation value. A direction of the plurality of directions having a maximum 
associated determined power spectral density value is then identified to thereby 
determine a direction of arrival for the ray. According to another trackmg aspect, a 
plurality of directions of arrival for the ray may be determined over a time interval 
20 such as a time slot or plurality of time slots. The determined plurality of directions of 
arrival may be processed with a tracking filter to produce a predicted direction of 

arrival for the ray. 

In other method aspects, a direction for an interfering radio signal is 
determined from the generated impairment autocorrelation matrix. A plurality of 

25 power spectral density values may be determined for the interfering radio signal at a 
like plurality of directions using the generated impairment autocorrelation matrix. A 
direction of the plurality of directions having a maximum associated determined 
power spectral density value may then be identified to thereby determine a direction 
for the interfering radio signal. The determined direction of arrival may be tracked 

30 over a time interval, such as a time slot or series of time slots. 

An apparatus for receiving a desired radio signal generated in a 
communications medium includes means, responsive to the commxmications medixmi, 
for receiving a plurality of radio signals from the communications medium. Means 
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are responsive to the means for receiving for generating a plurality of estimated signal 
medium response coefficients for a ray of the desired radio signal from the plurality of 
received radio signals. Means are also provided, responsive to the means for 
generating, for determining a direction for the desired radio signal from the plurality 
of estimated signal medium response coefficients for the ray of the desired radio 
signal. The means for determining a direction for the desired radio signal may 
comprise means for determining a direction of arrival for the ray from the 
corresponding plurality of estimated signal medium response coefficients for the ray, 
as well as means for determining a direction for the desired radio signal from the 
determined direction of arrival for the ray. 

The means for generating a plurality of estimated signal medium response 
coefficients may include means for generating a first plurality of estimated signal 
medium response coefficients. Means are responsive to the means for generating a 
first plurality of estimated signal medium response coefficients, for generating an 
impairment autocorrelation matrix from the plurality of received signals and the first 
plurality of estimated signal medium response coefficients. Additional means are 
provided, responsive to the means for generating an impairment autocorrelation 
matrix, for generating a branch metric from the generated impairment autocorrelation 
matrix and the first plurality of estimated signal medium response coefficients. 
Means are also provided, responsive to the means for generating a branch metric, for 
generating a symbol estimate from the generated branch metric. Means for generating 
a second plurality of estimated signal medium response coefficients from the plurality 
of received signals and the generated symbol estimate are provided. The apparatus 
may further include means for determining a direction for an interfering radio signal 
of the plurality of radio signals from the generated impairment autocorrelation matrix. 



Brief Description of the Drawings 

Some of the objects and advantages of the present invention having been 
stated, others will be more fully understood from the detailed description that follows 
and by reference to the accompanying drawings in which: 

FIG. 1 is a schematic illustration of an apparatus for determining a direction 
for a desired signal according to an embodiment of the present invention; 
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FIG. 2 is a schematic illustration of an apparatus for determining signal 
medium response coefficients and branch metrics using interference rejection 
combining according to an embodiment of the present invention; 

FIG. 3 is a schematic illustration of an apparatus for determining a direction 
5 for a desired signal according to an embodiment of the present invention; 

FIG. 4 is a schematic illustration of an apparatus for determining a direction 
for an interfering signal according to an aspect of the present invention; and 

FIG. 5 is a flov^chart illustration of operations for determining a direction for a 
desired signal according to an aspect of the present invention. 

Detailed Description of Preferred Embodiments 

The present invention now will be described more fully hereinafter with 
reference to the accompanying drawdngs, in which embodiments of the invention are 
shown. Those skilled in the art will appreciate that the invention may be embodied in 
many different forms and should not be construed as limited to the embodiments set 
forth herein; rather, these embodiments are provided so that this disclosure v^li be 
thorough and complete, and will fully convey the scope of the invention to those 
skilled in the art. In the drawings, like numbers refer to like elements throughout. 

Interference Rejection Combining 

A transmitted symbol sequence may be represented as s={s.l»".5Sl}. A signal 
r(t) presented at Na antennas may be take the vector form 

L(t) = £o(t; So) + ML(t), 
where a desired signal So(t;so) may be represented as 

L 

So(t;sJ- J^So^Mt^n^), 

w(t) is an impairment process, h(t,/T) is a composite signal medium, i.e., channel, 
response for the desired signal and T is the symbol interval. It may be assumed that 
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the impairment process contains both Ni undesired, e.g., interfering, signals and a 
zero-mean, white noise process n(t) such that 



w(/) = i;^/(^;^/)+2(0, (3) 

where the undesired signals Si(t;Si) have a similar form to the desired signal of 
equation (2) and wherein 

£{Z2(0«''(^)} = ^5(/-5X (4) 



10 where No is the noise power spectral density. A composite signal medium response 
for the desired signal may be resolved into two components: (1) a propagation signal 
medium response g(t, /T) and (2) a fixed transmit fiher response y(t, /T). The 
composite signal medium response h(t,/T) may then be assumed to have the form 

A— CO 

1 5 where t is time, Ts is the sample interval, M is the number of samples per symbol 
interval t, and gk are the channel tap coefficients for the channel response, the ^h 
channel tap is assumed to be constant over the transmit pulse. Those skilled in the art 
will appreciate that the above relationship may be used when the sample interval Ts is 
a fraction of the symbol interval T, corresponding to M equal to an integer greater 

20 than one. 

An equalizer may compute a log-likelihood metric LL(s;r(t)) for a 
hypothesized symbol s given a received signal r(t) according to 



25 



I L 



L 



\2 . 



2 



(6) 



where 
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r,= T g" y (7) 

A— 00 

represents the overall channel matched response, 
5 y,^lj\t-pT,yAt)dt (8) 

represents the filter matched response, 

CO 00 

10 

represents an interaction term for symbols si and 5^ and 

P(W)-(/-«)A. = U\t-ik^lMmf(t-U^mM)T,)dt (10) 

1 5 represents the pulse autocorrelation function. 

For a symbol-spaced equalizer withNyquist pulse shapes, equations (8) and 

(9) become 

y, = lj'it-pT)rSt)dt (11) 

20 and 



□0 00 



^1^1 S 5w^JtJt+/£jt+(/-m)Jfe+/' ^^^^ 

km 



respectively. 

In performing interference rejection combining, the impairment spatial 
25 autocorrelation matrix is estimated and inserted into the metric equations in the 
following manner: 
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15 



CO ^ CO ^ 



^=-00 



and 



where 



(15) 



In other words, the impairment autocorrelation matrix R is used to generate branch 
metrics which are compensated for the presence of interfering signals. The above 
approach is discussed in the aforementioned United States Patent Application Serial 
Number 08/284,775. 



Direction Estimation Using Phased Arrays 

Consider a received signal taking the form 

20 ^(0 = £(^>(0 + EiO. (1 6) 

where a(0) is the directional component of the signal presented at the array elements, 
s(t) is the transmitted signal and n(t) is a zero-mean, white noise component. Let 
y(t)=A"x(t) be the beamformed signal where A is the beamforming matrix. The 
25 spatial power spectrum of the signal x(t) may be estimated by computing 
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= (17) 

km] 



where the autocorrelation matrix Rxx is estimated by 



The spatial power spectrum Py(<l>) may similarly be computed as 



10 



where the autocorrelation matrix Ryy is estimated by 



R^(o)=^|;>:(/-«r,)/(/-«r,) 

= A^R„A- 



(20) 



1 5 The estimated direction is may then be computed by finding the value of the direction 
variable O that maximizes either equation (17) or (19). 

In the aforementioned United States Patent Application Serial Number 
08/655,930, a direction of arrival for a desired signal is estimated by maximizing 
power spectral density as shovra above using a interferer-compensated overall signal 

20 autocorrelation Rxx computed as Rxx = Rrr - Rii, where Rrr represents a running 

average of a received signal and Rii is an overall impairment autocorrelation matrix. 
Although this approach can provide an accurate estimate of direction of arrival, it may 
introduce noise which can reduce accuracy and/or which may require increased 
sampling to obtain a desired level of accuracy. 

25 The present invention arises from the realization that an accurate estimate of a 

direction for a desired signal may be obtained directly from estimated signal medium 
response coefficients (channel taps) for a ray of the desired signal, i.e., for a 
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component of the desired signal at one or more delays, instead of computing an 
overall autocorrelation for the desired signal and using it to determine a direction of 
arrival for the desired signal. According to the present invention, the estimated signal 
medium response coefficients gk,k+iM for a kth ray of the desired signal may be 
5 assumed to be a function of some direction parameter Oi<; for example, Ok=[<I'k,az,<l>k,ci] 
may describe the azimuth and elevation angles, respectively, of the medium response 
for the kth ray of the desired signal, where either 4>k,az or ^^,^1 may itself be known or 
assumed known. 

According to an aspect of the present invention, an autocorrelation for the kth 
1 0 ray may be computed as 



15 



allowing a spectral power density for the ^h ray to be computed as 

U''a''(0.)h^[A''£(S.)l 
|A"2(*.)|' 



The estimated direction for the kth ray can be estimated by finding the value of 
Ok that maximizes Pgk(^k). The gk,ic+iM represent a signal medium response for a kth 

20 ray, and in order to compute estimates of these components, the filtered data yk+iM 

may be used. In addition to this power spectral density maximization technique, other 
techniques may be used to determined direction for a ray based on the gk,k+iM, as will 
be discussed in greater detail below. 

Estimation of the direction for desired signal may computed from the 

25 determined direction of arrival of a ray of the desired signal in a number of ways. A 
particular ray might be selected for DOA computation because it may provide the 
greatest accuracy, for example. In other embodiments, a weighted combination of 
direction of arrival estimates for a plurality of rays of the desired signal may be 
computed to provide an accurate estimate of the direction of the source. 
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The direction of arrival information may be used for beamsteering and other 
applications. In addition, estimates of the directional components of an vmdesired 
signal may be generated using Ri,k+iM in a spatial power spectral density equation, and 
determining the direction of arrival for the undesired signal which maximizes power 
5 spectral density, as described above. This information may be advantageously used 
for antenna pattern nulling and other purposes, as described in greater detail below. 

Estimating Direction from Estimated Signal Medium Response Coefficients of a 
Component Ray 

10 

FIG. 1 illustrates the use of interference rejection combining (IRC) to produce 
estimated signal medium response coefficients g^ p according to an embodiment of the 
present invention. N antennas 110-1, , 110-N receive a plurality of radio signals 
from a composite signal 105 arriving at a receiving station 100. The plurality of 
15 received radio signals are passed to N radio units 115-1, - • • , 115-N. These signals 
are output from the radio units 115-1, • • 115-N and digitized by analog-to-digital 
(A/D) converters 120-1, • * * , 120-N. A synchronizer 125 performs synchronization 
on these signals, resulting in a N-dimensional vector output signal y^^ . Those skilled 

in the art will appreciate that synchronization may be performed jointly using all 
20 digitized signals to come up with one common sample timing, or may be performed 
individually for each digitized signal 

The signal y is passed to a branch metric determining means 130, which 

computes branch metrics n and hi,m for use in the symbol estimating means 150. The 
branch metric determining means 130 also produces estimated signal medium 
25 response coefficients g^^ p based on the received signal data 3^^ , as well as an 

impairment autocorrelation matrix Ri,p which is also estimated from the received data 
y . The symbol estimating means 150 produces an estimated symbol 5^, performing 

symbol hypothesis and detection using, for example, a Viterbi algorithm employed as 
part of a maximtmi likelihood sequence estimation (MLSE) process. A direction 
30 determining means 135 determines a direction for the desired signal from the 
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estimated signal medixim response coefficients gK.pfor the Ath ray of the desired 
signal. 

Those skilled in the art will appreciate that the radio units 115-1, , 115-N, 
A/D converters 120-1, • • • , 120-N, the synchronizer 125, the branch metric 

5 determining means 130, the symbol estimating means 150, and the direction 

determining means 135 may include a variety of commonly-used communications 
components such as filters, demodulators, digital signal processors and the like, the 
operation of which is well-know to those of skill in the art and need not be discussed 
in detail herein. In general, it will be appreciated that the functions of these 

10 components may be implemented using special-purpose hardware, software or 

firmware running on general or special purpose data processors, and combinations 
thereof. 

FIG. 2 illustrates an embodiment of the branch metric determining means 130 
of FIG. 1 in detail. A memory 220 stores estimated signal medium response 

1 5 coefficients g r-'i^gj^ ^ generated by signal medium response estimating means 
210. A delay 215 is provided between the signal medixmi response estimating means 
210 and the memory 220. The signal .medium response estimating means 210 
generates estimated signal medium response coefficients g^ ^ , • • • , i.N-hp ^^^^^ 
previously generated estimated signal medium response coefficients 

20 received fi:om the memory 220, received data and a symbol hypothesis , which 

may be a symbol estimate obtained using an MLSE or other symbol estimating 
process such as one performed by the symbol estimating means 150 depicted in FIG. 
1. To provide improved accuracy, a known training sequence, i.e., a synchronization 
sequence, may be employed to assist in estimation of the signal medium response in 

25 the signal medium response estimating means 210. For example, the training 

sequence may be used as the source of an initial symbol hypothesise, for determining 
the estimated signal medium response coefficients l;^^."^!;;^.,^^ • 

Signal medium response predicting means 240 uses the estimated signal 
medium response coefficients ^y-"^S^.^^^ ^o produce estimated fixture signal 

30 response medium coefficients g^^r- ^i^^., , > which are in turn used to compute the 
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branch metrics h/,m and r,. Impairment autocorrelation estimating means 225 
computes impairment autocorrelation values Ri.p (an N x N matrix) from an error 
signal Bp representing a difference between the received data and the symbol 
hypothesis s, . The impairment autocorrelation matrix Ri,p may be computed for 
5 sample time p and different delay indices i; setting i e {0} gives conventional 

interference rejection combining. A second delay 230 produces delayed impairment 
autocorrelation values Ri,p.i, for use in estimating the impairment autocorrelation 
values Ri,p. 

First branch metric computing means 245 computes a first branch metric hi,m , 
10 e.g., using equations (14) and (15), augmented by pulse autocorrelation values 

Pk-jHi-m)M provided by a pulse autocorrelation determining means 255, e.g., a look- 
up table of pulse autocorrelation values. A second branch metric computing 
determining means 250 computes a second metric n as in equation (13). It may be 
assumed that the summations in equation (13) are finite, i.e.. 



15 



N-l /-I 



*=0 



and that all terms are available. As illustrated in FIG. 1, the first and second metrics 
h|.n, and ri may be used by the symbol estimating means 150 to produce the estimated 
20 symbol values s, . Those skilled ui the art will appreciate that various apparatus may 
be used to produce the symbol values , such as MLSE sequence estimators, 
decision decoders and the like. 

The estimated signal mediimi response coefficients £o,p'"*'£a^-i,;'^°^ * 
desired signal and the estimated impairment autocorrelation values Ri,p may be used 

25 to determine a direction for a desired signal. FIG. 3 illustrates an exemplary structure 
in which a ray autocorrelation determining means 310 determines an autocorrelation 
Rgk,p for the kth ray of the desired signal across antennas, e.g., as described in 
equation (21). This value is used by a power spectral density determining means 320 
to determine spatial power spectral density values Pgk(<l>k) for the Alh ray as described 

30 in equation (22). Note that this may be performed for <l>k belonging to a set of 
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hypothesized directions {(X>k}. A ray direction determining means 330 determines a 
value of Ok that maximizes Pgk(C'k), producing an estimated direction of arrival 

6^^ for the kth ray. 

Those skilled in the art will appreciate that Pgk(<E>k) can be smoothed or 

5 averaged over time to give better estimates. Pgk(<I^k) can also be averaged over 
multiple antenna arrays to obtain spatial averaging. The estimated direction of a 
arrival 6^ for a kth ray, e.g., a main ray, may be used alone to estimate the direction 
of the desired signal, or a weighted combination of estimated directions of a arrival 
for a plurality of rays may be used to determine a direction for the desired signal. 

1 0 A direction of arrival for an interfering signal can also be obtained using the 

interference rejection combining process. Because the interference rejection 
combining process illustrated in FIG. 2 can provide impairment autocorrelation values 
Ri^p directly, the direction of arrival for the interfering signal may be directly 
determined by determining the interfering signal's power spectral density Pi(Ok) for a 

15 plurality of directions Ok in a power spectral density determining means 410, and 
detemiining the value 6,- that maximizes the interfering signal's power spectral 
density Pi(Ok) in an interfering signal direction determining means 420. 

Additional operations may also be performed to generate improve direction 
estimates. In a slotted time division multiple access (TDMA) system, for example, 

20 the du-ection estimates 6^^ and/or 6,. may be tracked over one or more slots using a 

tracking filter 350, 420 to yield predicted directions O^^ , 0, , Similarly, the 
estimated signal medium response coefficients go,p'"*'ljv-i,p impairment 
autocorrelation values Ri,p may be tracked within a slot or over a plurality of slots to 
produce improved direction estimates. 

25 Those skilled in the art will appreciate that the apparatus illustrated in FIGs. 2- 

4 may include a variety of commonly-used communications components. For 
example, the apparatus of FIG. 2 may be implemented using one or more digital 
signal processor (DSP) chips and/or application specific integrated circuits (ASICs). 
In general, it will be appreciated that the apparatus of FIGs. 2-4 may be implemented 

30 using special-purpose hardware, software or firmware running on general or special 
purpose data processors, and combinations thereof. 
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It will be appreciated that techniques for determining directions for desired 
and interfering signals other than the aforementioned spectral power density 
maximization approach may be used with the present invention. For example, a 
technique could be utilized which determine direction using signal strength estimates 
5 computed from the gk , similar to the technique described in United States Patent 
Application Serial Number 08/597,073. Other approaches, such ones based on the 
MUSIC techniques described in "A Signal Subspace Approach to Multiple Emitter 
Location and Spectral Estimation", by Schmidt (Stanford University Ph.D. 
dissertation, 1981), or the ESPRIT techniques described in "ESPRIT, Estimation of 
1 0 Signal Parameters via Rotational Invariance Techniques", by Roy (Stanford 
University Ph.D. dissertation, 1987), could be used as well. 

The directional information for desired and interfering signals may be used for 
a number of different purposes. For example, the estimated directions and O. 
may be used to steer maximum gain portions of an antenna pattern to advantageously 
15 receive a desired signal, in a manner similar to the beamsteering operations described 
in the aforementioned U.S. Patent Application Serial No. 08/655,930. Estimates of 
interfering signal directions may be used, for example, for nulling a receive antenna 
pattem to reduce interference, as well as for such features as predicting mobile unit 
handoffs and the like. 

20 FIG. 5 is a flowchart illustration of operations for determining a direction for a 

desired signal arriving at a receiving station (Block 500). A plurality of radio signals 
are received at a plurality of antennas (Block 510). A plurality of estimated signal 
medium response coefficients for a ray of the desired radio signal are generated from 
the plurality of received signals (Block 520), e.g., using in interference rejection 

25 combining technique as described with respect to FIG. 2. A direction for the desired 
signal is determined from the plurality of estimated signal medium response 
coefficients for the ray of the desired signal (Block 530), for example, by computing 
an autocorrelation for the ray and determining a direction that maximizes a power 
spectral density function of the autocorrelation. 

30 In the drawings and specification, there have been disclosed typical 

embodiments of the invention and, although specific terms are employed, they are 
used in a generic and descriptive sense only and not for purposes of limitation, the 
scope of the invention being set forth in the following claims. 



wo 99/33197 



17 



PCT/US98/27104 



THAT WHICH IS CLAIMED IS: 

1 . A method of determining a direction for a desired radio signal 
generated in a commimications medium, the method comprising the steps of: 

receiving a plurality of radio signals from the communications medium; 
generating a plurality of estimated signal medium response coefficients for a 
5 ray of the desired radio signal from the plurality of received radio signals; and 
determining a direction for the desired radio signal from the plurality of 
estimated signal medium response coefficients for the ray of the desired radio signal. 

2. A method according to Claim 1 , wherein said step of receiving 
comprises the step of receiving the plurality of radio signals at a like plurality of 
antennas. 

3. A method according to Claim 2: 

wherein said step of determining a direction for the desired radio signal 
comprises the step of determining a direction of arrival for the ray from the 
corresponding pluraUty of estimated signal medium response coefficients for the ray; 
5 and 

wherein said step of determining a direction for the desired radio signal 
comprises the step of determining a direction for the desired radio signal from the 
determined direction of arrival for the ray. 

4. A method according to Claim 2: 

wherein said step of generating a plurality of estimated signal medium 
response coefficients comprises the steps of: 

generating a series of pluralities of estimated signal medium response 
5 coefficients over a time interval; and 

processing the series of pluralities of estimated signal medium 
response coefficients with a tracking filter to produce a plurality of predicted 
signal medium response coefficients; and 

wherein said step of determining a direction for the desired radio signal 
10 comprises the step of determining a direction for the desired radio signal from the 



wo 99/33197 PCT/US98/27104 

18 

plurality of predicted signal medium response coefficients for the ray of the desired 
radio signal. 

5. A method according to Claim 4, wherein the desired radio signal is 
formatted according to a slotted commimications protocol defined according to a 
plurality of time Slots, wherein said step of generating a series of pluralities of 
estimated signal medium response coefficients over a time interval comprises the step 

5 of generating a series of pluralities of estimated signal medium response coefficients 
over at least one of said plurality of time slots. 

6. A method according to Claim 2, wherein the desired radio signal 
represents a sequence including a predetermined training sequence, and wherein said 
step of generating a plurality of estimated signal medium response coefficients 
comprises the step of generating a plurality of estimated sigiial medium response 

5 coefficients from the plurality of received radio signals based on the predetermined 
training sequence. 

7. A method according to Claim 2, wherein said step of generating a 
plurality of estimated signal medium response coefficients for the ray comprises the 
steps of: 

generating a first plurality of estimated signal medium response coefficients; 
5 generating an impairment autocorrelation matrix from the plurality of received 

signals and a symbol estimate; 

generating a branch metric from the generated impairment autocorrelation 
matrix and the first plurality of estimated signal medium response coefficients; 
generating a new symbol estimate from the generated branch metric; and 
10 generating a second plurality of estimated signal medium response coefficients 

from the plurality of received signals and the generated new symbol estimate. 



8. A method according to Claim 7: 

wherein said step of generating an impairment autocorrelation matrix 
comprises the steps of: 
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generating a plurality of impairment autocorrelation matrices over a 

5 time interval; and 

processing the generated plurality of impairment autocorrelation 
matrices with a tracking filter to produce a predicted impairment 
autocorrelation matrix; and 

wherein said step of generating a branch metric comprises the step of 
1 0 generating the branch metric from the predicted impairment autocorrelation matrix. 

9. A method according to Claim 8, wherein the desired radio signal is 
formatted according to a slotted communications protocol defined according to a 
plurality of time slots, and wherein said step of generating a plurality of impairment 
autocorrelation matrices over a time interval comprises the step of generating a 

5 plurality of impairment autocorrelation matrices over at least one time slot of said 
plurality of time slots. 

10. A method according to Claim 7, wherein said step of generating a 
symbol estimate comprises the step of generating a symbol estimate from the 
generated branch metric according to a maximum likelihood sequence estimation 
(MLSE) process. 

11. A method according to Claim 3 , wherein said step of determining a 
direction of arrival for the ray comprises the steps of: 

determining an autocorrelation value for the ray of the desired radio signal; 

determining a plurality of power spectral density values for the ray at a like 
5 plurality of directions using the determined autocorrelation value; and 

identifying a direction of the plurality of directions having a maximum 
associated determined power spectral density value to thereby determine a direction 
of arrival for the ray. 



12. A method according to Claim 3 : 

wherein said step of determining a direction of arrival for the ray comprises 
the steps of: 
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determining a plurality of directions of arrival for the ray over a time 

5 interval; and 

processing the determined plurality of directions of arrival with a 
tracking filter to produce a predicted direction of arrival for the ray; and 
wherein said step of determining a direction for the desired radio signal from 
the determined direction of arrival for the ray comprises the step of determining a 
1 0 direction for the desired radio signal from the predicted direction of arrival for the ray. 

13. A method according to Claim 12, wherein the desired radio signal is 
formatted according to a slotted communications protocol defined according to a 
plurality of time slots, and wherein said step of determining a plurality of directions of 
arrival for the ray over a time interval comprises the step of determining a plurality of 

5 directions of arrival for the ray over at least one of said plurality of time slots. 

14. A method according to Claim 7, further comprising the step of 
determining a direction for an interfering radio signal of the plurality of radio signals 
from the generated impairment autocorrelation matrix. 

15. A method according to Claim 8, further comprising the step of 
determining a direction for an interfering radio signal of the plurality of radio signals 
from the predicted impairment autocorrelation matrix. 

16. A method according to Claim 9, further comprising the step of 
determining a direction for an interfering radio signal of the plurality of radio signals 
from the predicted impairment autocorrelation matrix. 

17. A method according to Claim 14, wherein said step of determining a 
direction for an interfering radio signal comprises the steps of: 

determining a plurality of power spectral density values for the interfering 
radio signal at a like plurality of directions using the generated impairment 
5 autocorrelation matrix; and 
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identifying a direction of the plurality of directions having a maximum 
associated determined power spectral density value to thereby determine a direction 
for the interfering radio signal. 



1 8. A method according to Claim 17, wherein said step of determining a 
direction for the interfering radio signal comprises the steps of: 

determining a plurality of directions of arrival for the interfering radio 
signal over a time interval; and 
5 processing the determined plurality of directions of arrival with a 

tracking filter to produce a predicted direction of arrival for the interfering 
radio signal. 

19. A method according to Claim 1 8, wherein the desired radio signal is 
formatted according to a slotted conununications protocol defined according to a 
plurality of time slots, and wherein said step of determining a plurality of directions of 
arrival for the interfering radio signal over a time interval comprises the step of 

5 determining a plurality of directions of arrival for the interfering radio signal over at 
least one of said time slots. 

20. An apparatus for receiving a desired radio signal generated in a 
communications medium, the station comprising: 

means, responsive to the communications medium, for receiving a plurality of 
radio signals from the communications medium; 
5 means, responsive to said means for receiving, for generating a plurality of 

estimated signal medium response coefficients for a ray of the desired radio signal 
from the plurality of received radio signals; and 

means, responsive to said means for generating, for determining a direction for 
the desired radio signal from the plurality of estimated signal medium response 
1 0 coefficients for the ray of the desired radio signal. 

21 . An apparatus according to Claim 20, wherein said means for receiving 
comprises a plurality of antennas, a respective one of which produces a respective one 
of the plurality of received signals. 
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22. An apparatus according to Claim 21 : 

wherein said means for determining a direction for the desired radio signal 
comprises means for determining a direction of arrival for the ray from the 
corresponding plurality of estimated signal medium response coefficients for the ray; 
5 and 

wherein said means for determining a direction for the desired radio signal 
comprises means for determining a direction for the desired radio signal from the 
determined direction of arrival for the ray. 

23 . An apparatus according to Claim 2 1 : 

wherein said means for generating a plurality of estimated signal medium 
response coefficients comprises: 

means for generating a series of pluralities of estimated signal medium 
5 response coefficients over a time interval; and 

means for processing the series of pluralities of estimated signal 
medium response coefficients with a tracking filter to produce a plurality of 
predicted signal medium response coefficients; and 

wherein said means for determining a direction for the desired radio signal 
10 comprises means for determining a determining a direction for the desired radio signal 
from the plurality of predicted* signal medium response coefficients for the ray of the 
desired radio signal. 

24. An apparatus according to Claim 23, wherein the desired radio signal 
is formatted according to a slotted communications protocol defined according to a 
plurality of time slots, wherein said means for generating a series of pluralities of 
estimated signal medium response coefficients over a time interval comprises means 

5 for generating a series of pluralities of estimated signal medium response coefficients 
over at least one of said plurality of time slots. 

25. An apparatus according to Claim 22, wherein the desired radio signal 
represents a sequence including a predetermined training sequence, and wherein said 
means for generating a plurality of estimated signal medium response coefficients 
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comprises means for generating a plurality of estimated signal medium response 
5 coefficients from the plurality of received radio signals based on the predetermined 
training sequence. 

26. An apparatus according to Claim 23, wherein said means for 
generating a plurality of estimated signal medium response coefficients for the ray 
comprises: 

means for generating a first plurality of estimated signal medium response 
5 coefficients; 

means, responsive to said means for generating a first plurality of estimated 
signal medium response coefficients, generating an impairment autocorrelation matrix 
from the plurality of received signals and a symbol estimate; 

means, responsive to said means for generating an impairment autocorrelation 
1 0 matrix, for generating a branch metric from the generated impairment autocorrelation 
matrix and the first plurality of estimated signal medium response coefficients; 

means, responsive to said means for generating a branch metric, for generating 
a new symbol estimate from the generated branch metric; and 

means, responsive to said means for generating a symbol estimate, for 
1 5 generating a second plurality of estimated signal medium response coefficients from 
the plurality of received signals and the generated new symbol estimate. 

27. An apparatus according to Claim 26: 

wherein said means for generating an impairment autocorrelation matrix 
comprises: 

means for generating a plurality of impairment autocorrelation 
5 matrices over a time interval; and 

means, responsive to said means for generating a plurality of 
impairment autocorrelation matrices, for processing the generated plurality of 
impairment autocorrelation matrices with a tracking filter to produce a 
predicted impairment autocorrelation matrix; and 
10 wherein said means for generating a branch metric comprises means for 

generating the branch metric from the predicted impairment autocorrelation matrix. 
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28. An apparatus according to Claim 27, wherein the desired radio signal 
is formatted according to a slotted communications protocol defined according to a 
plurality of time slots, and wherein said means for generating a plurality of 
impairment autocorrelation matrices over a time interval comprises means for 

5 generating a plurality of impairment autocorrelation matrices over at least one time 
slot of said plurality of time slots. 

29. An apparatus according to Claim 26, wherein said means for 
generating a symbol estimate comprises means for generating a symbol estimate from 
the generated branch metric according to a maximum likelihood sequence estimation 
(MLSE) process. 

5 

30. An apparatus according to Claim 22, wherein said means for 
determining a direction of arrival for the ray comprises: 

means for determining an autocorrelation value for the ray of the desired radio 

signal; 

10 means, responsive to said meaiis for determming an autocorrelation value for 

the ray, for determining a plurality of power spectral density values for the ray at a 
like plurality of directions using the determined autocorrelation value; and 

means, responsive to said means for determining a plurality of power spectral 
density values, for identifying a direction of the plurality of directions having a 

1 5 maximum associated determined power spectral density value to thereby determine a 
direction of arrival for the ray. 

31. An apparatus according to Claim 22: 

wherein said means for determining a direction of arrival for the ray 
comprises: 

means for determining a plurality of directions of arrival for the ray 

5 over a time interval; and 

means, responsive to said means for determining a plurality of 
directions of arrival, for processing the determined plurality of directions of 
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arrival with a tracking filter to produce a predicted direction of arrival for the 
ray; and 

1 0 wherein said means for determining a direction for the desired radio signal 

from the determined direction of arrival for the ray comprises means for determining 
a direction for the desired radio signal from the predicted direction of arrival for the 
ray. 

32. An apparatus according to Claim 3 1 , wherein the desired radio signal 
is formatted according to a slotted communications protocol defined according to a 
plurality of time slots, and wherein said means for determining a plurality of 
directions of arrival for the ray over a time interval comprises means for determining 

5 a plurality of directions of arrival for the ray over at least one of said plurality of time 
slots. 

33. An apparatus according to Claim 27, fiirther comprising means for 
determining a direction for an interfering radio signal of the plurality of radio signals 
from the generated impairment autocorrelation matrix. 

34. An apparatus according to Claim 27, further comprising means for 
determining a direction for an interfering radio signal of the plurality of radio signals 
from the predicted impairment autocorrelation matrix. 

35. An apparatus according to Claim 28, further comprising means for 
determining a direction for an interfering radio signal of the plurality of radio signals 
from the predicted impairment autocorrelation matrix. 

36. An apparatus according to Claim 33, wherein said means for 
determining a direction for an interfering radio signal comprises: 

means for determining a plurality of power spectral density values for the 
interfering radio signal at a like plurality of directions using the generated impairment 
5 autocorrelation matrix; and 

means, responsive to said means for determining a plurality of power spectral 
density values, for identifying a direction of the plurality of directions having a 
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maximum associated deteraiined power spectral density value to thereby determine a 
direction for the interfering radio signal. 



37. An apparatus according to Claim 36, wherein said means for 
determining a direction for the interfering radio signal comprises: 

means for determining a plurality of directions of arrival for the 
interfering radio signal over a time interval; and 
5 means, responsive to said means for determining a plurality of 

directions of arrival, for processing the determined plurality of directions of 
arrival with a tracking filter to produce a predicted direction of arrival for the 
interfering radio signal. 

38. An apparatus according to Claim 37, wherein the desired radio signal 
is formatted according to a slotted communications protocol defined according to a 
plurality of time slots, and wherein said means for determining a plurality of 
directions of arrival for the interfering radio signal over a time interval comprises 

5 means for determining a plurality of directions of arrival for the interfering radio 
signal over at least one of said time slots. 
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